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The dynamical behavior of a weakly damped harmonic 
chain in a spatially periodic potential (Frenkel-Kontorova 
model) under the subject of an external force is investigated. 
We show that the chain can be in a spatio-temporally chaotic 
state called fluid-sliding state. This is proven by calculating 
correlation functions and Lyapunov spectra. An effective tem- 
perature is attributed to the fluid-sliding state. Even though 
the velocity fluctuations are Gaussian distributed, the fluid- 
sliding state is clearly not in equilibrium because the equipar- 
tition theorem is violated. We also study the transition be- 
tween frozen states (stationary solutions) and molten states 
(fluid-sliding states). The transition is similar to a first-order 
phase transition, and it shows hysteresis. The depinning- 
pinning transition (freezing) is a nucleation process. The 
frozen state contains usually two domains of different par- 
ticle densities. The pinning-depinning transition (melting) is 
caused by saddle-node bifurcations of the stationary states. 
It depends on the history. Melting is accompanied by precur- 
sors, called micro-slips, which reconfigurate the chain locally. 
Even though we investigate the dynamics at zero tempera- 
ture, the behavior of the Frenkel-Kontorova model is quali- 
tatively similar to the behavior of similar models at nonzero 
temperature. 

PACS numbers: 46.10.+Z, 46.30.Pa, 68.35.Rh 



I. INTRODUCTION 

Systems with many degrees of freedom which are 
pinned in some external potential are very common in 
condensed matter. Examples are fluid-fluid interfaces in 
porous media QJ^] , flux-lattices in type-II superconduc- 
tors ||, and charge density waves B to mention only a 
few. Also dry friction (i.e., solid-solid friction) belongs to 
this class of systems because the asperities of the surfaces 
interlock. 

A common feature of all these systems is a strongly 
nonlinear mobility. If one applies some field or force F on 
the system, the mobility is zero below some usually well- 
defined threshold F c . Above this threshold the mobility is 
nonzero. In general it is some nonlinear function of the 
applied force F. The transition from a pinned system 
with zero mobility to a depinned one with some finite, 
nonzero mobility is called the pinning-depinning transi- 
tion. It can be understood as a kind of "melting" which 
happens far from thermal equilibrium. The process is a 
typical nonequilibrium one because of two reasons. First, 
there is no ground state for F ^ and the pinned sys- 
tem has to be in some metastable state. Due to thermal 
fluctuations the system can overcome the barrier of the 
metastable state and move into another metastable state 
with less energy. This phenomenon leads to creeping 



with a very low mobility. Second, beyond the pinning- 
depinning transition energy flows through the system at 
a constant rate which is given by the mobility times F 2 . 
This flow is usually not small. Thus it cannot be de- 
duced from linear response theory which works only near 
thermal equilibrium. The mobility of the sliding state 
strongly depends on the kind of energy dissipation. 

The inverse process of this nonequilibrium melting 
is the depinning-pinning transition which is a kind of 
nonequilibrium "freezing" . Both kind of transitions do 
not have to occur at the same value of the applied force 
F. The behavior depends strongly on whether the de- 
grees of freedom (i.e., flux lines, atoms etc.) have inertia 
or inertia is negligible compared to dissipative forces. In 
the case of strong dissipation the motion is overdamped. 
The pinning-depinning transition is in most cases of 
second order and indistinguishable from the depinning- 
pinning transition. Typical examples for such systems 
are flux lines in type-II superconductors and charge- 
density waves. If the motion is underdamped, hysteresis 
is possible because the inertia can overcome a pinning 
center. This is intuitively clear if one imagines the sim- 
plest model system of this kind, namely a particle in a 
spatially periodic potential || . 

Another important aspect of the collective behavior 
of pinned systems is whether the potential caused by the 
pinning centers is regular or irregular (quenched random- 
ness). Often the pinning landscape is random. This case 
together with a purely dissipative diffusion- like dynamics 
has been studied extensively in the literature j6jQ]. 

The aim of this paper is to study the opposite case 
in a fairly simple model, namely, the Frenkel-Kontorova 
(FK) model |8|. There is no quenched randomness. All 
pinning centers are identical forming a regular array. Fur- 
thermore, all pinned objects are identical and have a 
mass. The damping is assumed to be weak. We will see 
that weak damping is responsible for randomness that is 
caused by chaotic motion. Important physical applica- 
tions of the FK model are arrays of identical Josephson 
junctions S and adsorbate layers on clean crystal sur- 
faces @. 

In this paper we consider the one-dimensional FK 
model. The equation of motion (in dimensionless units) 
reads 



JXj 



j + i — 2xj — bsiriXj + F, (1) 



where Xj is the position of particle j, 7 is the damping 
constant, b is the strength of the external potential, and 
F is the external force. The time derivative is denoted 
by a dot. In order to avoid effects due to the boundary 
layers we choose periodic boundary conditions, i.e., 



Xj+ N = Xj + 2irM, 



(2) 



1 



where N is the number of particles and M is an arbitrary 
integer. The periodic boundary conditions fixes the av- 
erage particle distance a to a = 2ttM/N. Because of 
symmetry a can be restricted to interval [0, tt] without 
loss of generality.. 

Together with the previous paper |p"l|| in which we have 
already investigated periodic and quasiperiodic solutions, 
the aim is to give a detailed investigation of the dynam- 
ical behavior in the weakly damped case for long chains 
(i.e., N > 100). 
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FIG. 1. The velocity-force characteristic and the effective 
temperature of the fluid-sliding state. The different branches 
belong to stationary states (SS), fluid-sliding states (FSS), 
and solid-sliding states (SSS). In the simulations the applied 
force F was decreased (squares and dotted lines) or increased 
(triangles and solid lines) with a constant rate (|-F| = 10~ 7 ). 
The velocity at each data point is the average over a time 
interval of 10 4 time units. The upper inset shows 20 hysteresis 
loops between SS and FSS from a simulation where F was 
moved in the interval [0.09, 0.15] forward and backward at a 
rate of \F\ — 10~ 6 . The lower inset shows the particle velocity 
distribution for F = 0.26. The solid line shows the fit of the 
data point with a Gaussian. The parameters are N = 233, 
M = 89, b = 2, and 7 = 0.05. 

In this paper we deal with spatio-temporal chaotic so- 
lutions (called fluid-sliding states) and the transition be- 
tween these solutions and the stationary states. The typ- 
ical behavior is summarized in figure 0. Fig. [jja) shows 
the velocity-force characteristic. We see hysteresis loops 
between three different branches which belong to differ- 



ent types of solutions. The states with the largest average 
sliding velocities v are the solid-sliding states which are 
characterized by a chain with nearly no internal vibra- 
tions. These states have the maximum possible mobility, 
i.e., I/7. The solid-sliding state becomes unstable due to 
first-order parametric resonance if v is below some criti- 
cal value pJJJ , 

The second type of sliding state is the fluid-sliding 
state. In general one can summarize all sliding states 
having not the maximum mobility into this category. 
But strictly speaking the name makes only sense if these 
states are spatio-temporally chaotic. For larger values 
of the damping constant this is not the case as we have 
seen in the previous paper. The chaotic vibration in the 
fluid-sliding state can be characterized by an effective 
temperature. The temperature of the fluid-sliding states 
of Fig. 0(a) is shown in figure 0(b). Even though the dis- 
tribution of the particle velocity ij is gaussian [see inset 
of Fig. 0(b)] we will show that the fluid-sliding state is a 
nonequilibrium state. A very specific test to show this is 
the violation of the equipartition theorem for the phonon 



modes (see Sec. II A) 



The third type of states are the stationary ones. Their 
mobility is zero. In order to model the creeping due to 
thermal activation one has to add a white-noise term to 
the equation of motion (0). We have not done this be- 
cause the qualitative behavior does not change very much 
as long as the thermal energy is much smaller than the 
amplitude of the external potential. This is confirmed in 
numerical simulations of similar models po| , p^jl3|| . For 
example, the hysteresis seen in Fig. (a) still exist for 
nonzero but small temperatures 10 1% . For this behav- 
ior it seems to be important that the system has many 
degrees of freedom because in the case of N = 1 the 
hysteresis disappears even for infinitesimal small noise 
amplitude ||. 

Fig. 0(b) clearly shows that nonequilibrium melting 
and freezing is accompanied by an abrupt change of 
the temperature of the chain. The transition is like a 
first-order one in thermal equilibrium but the pinning- 
depinning transition point is larger than the depinning- 
pinning transition point. Thus hysteresis occurs. The 
transition points fluctuate [see inset of Fig. 0(a)], espe- 
cially the pinning-depinning transition point. 

The paper is organized as follows: In section |l| we 
investigate in detail the fluid-sliding state. We show 
that it is indeed spatio-temporally chaotic. For chains 
with a/2ir near an integer value we found a pronounced 
transition from a kink-dominated sliding state and the 
fluid-sliding state. This transition is relatively sharp even 
though there is no hysteresis. But it becomes hysteretic 
for small N. The depinning-pinning transition and the 



pinning-depinning transition are discussed in section III. 
We show that nonequilibrium freezing is similar to ordi- 
nary freezing whereas melting is clearly different. The 
pinning-depinning transition point depends on the sta- 
tionary state. Furthermore local rearrangements (micro- 
slips) of the chain may occur before the transition. In 
section 



models. 



IV we compare our results with results of similar 
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FIG. 2. An example for spatio-temporal chaos. Each solid 
line is a snapshot of the system. The time interval between 
two successive snapshots is St = 4n/v (v is the average sliding 
velocity). A particular snapshot is highlighted. The parame- 
ters are N = 233, M = 89, b = 2, 7 = 0.05, and F = 0.14. 



II. CHAOTIC SLIDING 



A. Spatio-temporal chaos 

Fig ^ is of course not a proof that the chain slides 
chaotically. It is well-known that chaotic motion is char- 
acterized by the sensitivity on the initial conditions. It is 
measured by the largest Lyapunov exponent A max , which 
is the rate of divergence (or convergence, if it is neg- 
ative) of trajectories in phase space that start out in- 
finitely close to each other J L| . Figure || shows that the 
fluid-sliding states in Figs. [2 and |l| arc indeed temporally 
chaotic. 



U 




Decreasing the damping constant 7 increases the com- 
plexity of the sliding state from periodic motion via 
quasi-periodic motion (which is usually spatially chaotic, 
see the preceding paper) to spatio-temporal chaos. An 
example of the latter one is shown in figure |2[ 

The aim of this section is to investigate and to char- 
acterize the chaotic-sliding state which we call the fluid- 
sliding state. First of all we see in Fig. |l|(a) that the 
velocity-force characteristic of this state is nearly struc- 
tureless. This has to be compared with the case of peri- 
odic and quasiperiodic motion where a multitude of hys- 
teresis loops appear (see the preceding paper). Here there 
occur only hysteresis loops between the solid-sliding state 
(where the particles are shaken so fast that they nearly 
do not feel the external potential), the fluid-sliding state 
and the stationary states. 
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FIG. 3. The maximum Lyapunov exponent of the 
fluid-sliding state as a function of the applied force F. The 
parameters are N = 144, M = 55, b = 2, and 7 = 0.05. 



FIG. 4. The normalized velocity correlation function Cj 
of the fluid-sliding state. To guide the eye the numerical re- 
sults (denoted by squares) are connected by a solid line. The 
inset shows a logarithmic plot of \cj\. The straight line is 
the function exp(— j/2). The parameters are the same as in 
Fig. | 

But is it also spatially chaotic? In order to answer 
this question we have calculated the normalized velocity 
correlation function Cj defined by 



where 



a, 



((/,)) 



_ {{XiX l+j ) 



to; 



to 



1 r i N 

lim - / — V fi(t)dt 

.7 = 1 



(3) 



(4) 



For the same parameters as in Fig. g the result is shown 
in figure ||. One clearly sees that Cj is a rapidly decay- 
ing oscillatory function. The envelope seems to be pro- 
portional to exp(— with a correlation length (k2. 
Because of N ^> £ and A max > the fluid-sliding state is 
spatio-temporally chaotic. 

A very strong criterion for spatio-temporal chaos is 
that the number of positive Lyapunov exponents is pro- 
portional to TV for large N. We have calculated the Lya- 
punov spectrum with the method described in Rcf. |lq | 
for various values of N. Figure ^ shows the cumulative 
density pjv(A) of Lyapunov exponents, i.e., the probabil- 
ity to find a Lyapunov exponent larger than A. The result 
is typical for spatio-temporal chaos [HJ . In the thermo- 
dynamic limit (i.e., N — > 00) the sequence of cumulative 
densities converges uniformly to p^. Thus, for large 
N the number of positive Lyapunov exponents is indeed 



3 



proportional to N. Fig. |5| shows clearly that the spatio- 
temporally chaotic nature of the fluid-sliding state does 
not depend on the commensurability of the chain. 




0.30 



FIG. 5. Lyapunov spectra for (a) the most commensurate 
case (i.e. a = 0) and (b) the most incommensurate case (i.e. 
a/2% -» (3- V5)/2 = 2- golden mean). The spectra for two 
different system sizes are shown. Squares and solid lines de- 
note (a) N = 20, (b) N = 13, and (a) N = 200, (b) N = 233. 
The other parameters are b — 2, 7 = 0.05, and (a) F — 0.13, 
(b)F = 0.3. 

Because the fluid-sliding state is spatio-temporally 
chaotic it makes sense to introduce an effective tempera- 
ture. In the dimensionless units of the equation of motion 
([j]) it is the average kinetic energy in the frame co-moving 
with the center of mass, i.e., 



T 



where v is the average sliding velocity 



(5) 



(6) 



In the previous paper we have derived a formula for the 
applied force F in terms of the first and second moments 
of the particle velocity [Eq. (6) in Ref. jjl)]. With the 
help of this formula the temperature can be expressed in 
terms of the applied force and the average sliding veloc- 
ity: 



(7) 



Figure 0(b) shows the temperatures of the fluid-sliding 
states of the velocity- force characteristic in Fig. |l|(a). 

Even though the temperature of the solid-sliding states 
and the stationary states is formally zero in accordance 
with (0), it does not make sense to call @ a "temper- 
ature" in regular, non-chaotic sliding states or station- 
ary states. The periodic and quasi-periodic domain-like 
states, for example, investigated in the previous paper 
have also non-zero "temperature" . 

In the frame co-moving with the center of mass, the 
chain is shaken by the washboard wave (i.e., the external 
potential) and moves in a spatio-temporally chaotic way. 
The gaussian distributed velocities might suggest that 
the chain is in thermal equilibrium. But is it true? This 
raises the following question of general interest: Can we 



replace the spatio-temporally chaotic chain by an equiv- 
alent system which is in thermal equilibrium? Or more 
general, is it possible to describe the chaotic attractor 
of a weakly damped and strongly driven Hamiltonian 
system with many degrees of freedom (infinitely many 
in the thermodynamic limit) by an equivalent undriven 
and undamped system? As a consequence of a positive 
answer one would expect that the equipartition theo- 
rem from thermodynamics holds, i.e. the ensemble av- 
erages of qjdH/dqj and PjdH/dpj are independent of 
j [H(qi,p±, . . . ,qj,pj, . . .) is the Hamilton function, the 
<7j's are the generalized coordinates, and the pj's are the 
corresponding canonical momenta]. In numerical simu- 
lations one usually replaces the ensemble average by the 
temporal average assuming that the ergodicity hypothe- 
sis holds. An obvious candidate for a test of the equipar- 
tition theorem would be the particle momentum in the 
co-moving frame, i.e., Xj—v. But this is not a wise choice: 
because of symmetry the result has to independent of j. 
A better choice is its spatial Fourier transform, i.e., 



1 N 



2iirkj/N 



k = 0,l,...,N -1. (8) 



That is, we want to check whether the average kinetic 
energy 



ejj = lim — 

r — >oo T 



\pk(t)\ 2 dt 



(9) 



of the phonon modes is equipartitioned or not. Figu re H 
shows that the equipartition theorem is not fulfilled p7| T 
This is a clear signature for the fact that the fluid- sliding 
state is a state far away from thermal equilibrium. There- 
fore it is not possible to develop a theory for this state 
based on equilibrium thermodynamic. The violation of 
the equipartition theorem is equivalent with non-zero ve- 
locity correlations for j ^ because e/- is the modulus of 
the Fourier transform of Ca . 
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FIG. 6. The average kinetic energy et of the phonon 
modes for three different values of the applied force F. The 
parameters are the same as in Fig. |lj 



B. The transition between fluid-sliding state and 
kink-dominated sliding state 

When a/2-K approaches an integer value the velocity- 
force characteristic of the fluid-sliding state develops a 
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relatively sharp transition step at a characteristic value 
of the applied force F. An example for a/2n = 1/20 
is shown in figure [?]. For long enough chains no hys- 
teresis is observable. For small chains we get bistability 
between different types of sliding states. Similar results 
has been found in a generalized FK model by Braun et al. 
pO| , ^4|Jl8| . The aim of this section is to answer the follow- 
ing obvious questions: What is the nature of the different 
sliding states? Why does the bistability depends on N? 
Can we understand this transition and where does it oc- 
cur? 
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FIG. 7. The transition between the fluid-sliding state and 
the kink-dominated sliding state. The velocity-force charac- 
teristics for N = 500 (filled symbols) and N = 200 (open 
symbols) are shown. Squares and dotted lines (triangles 
and solid lines) indicate decreasing (increasing) applied force 
F. The rate |F| is always 10~ 7 except for N — 500 and 
F G (0.26,0.29) where it is 10" 8 . The dashed line indicates 
the sound velocity which is equal to a. The parameters are 
a = tt/10, b = 2, and 7 = 0.05. 

First we take a more detailed look at the dynamics 
below and above the transition (see Fig. ||). The motion 
far below the threshold is almost regular. It corresponds 
to one of the multi-domain states we have discussed in the 
previous paper. There are two domain types, a stationary 
one with a = [it is responsible for the tilted lines in 
Fig. ||(a)] and a sliding one. Often the sliding domains 
are so small that they are actually 27r-kinks, and larger 
sliding domains can be interpreted as clusters of 27r-kinks 
p0| . That is, multi-domain states like the example of 
Fig. ||(a) are nonuniform distributions of 27r-kinks. Thus, 
we call this state the kink- dominated sliding state. Note, 
that there are M kinks but no antikink. 

It is well-known that kinks (and antikinks) cannot 
travel faster than the sound velocity (which is equal to 
one in our case). Each kink or antikink needs therefore 
at least N time units to travel through the whole chain. 
After that time a chain with M kinks will be shifted by 
2irM. Therefore, the average sliding velocity of a state 
like the one of Fig. ||(a) has to be less than 2ttM/N = a. 



Fig. [7] shows that it is actually much below the sound 
velocity. 
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FIG. 8. The dynamics of the kink-dominated sliding state 

(a) and the fluid-sliding state (b). Several snapshots are 
shown taken at equidistant time steps (a) St — 2n/v and 

(b) St — 20-k/v. In each case a particular snapshot is high- 
lighted. The parameters are (a) F = 0.2 and (b) F = 0.3 and 
N = 500, M = 25, (i.e., a = tt/10), b = 2, and 7 = 0.05. 
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FIG. 9. The average sliding velocity v (averaged over 
intervals of 200 time units) as function of time. The velocity 
vo of the quiescent phase is denoted by a dotted line. The 
parameters are N = 500, M = 25, (i.e., a — tt/10), b = 2, 
7 = 0.05, and F = 0.24. 

If the transition point is approached from below, the 
behavior depends on whether the chain is long or short. 
For long chains with many kinks, the average sliding ve- 
locity v starts to increase with F faster and faster. Later 



5 



on the increase slows down. We define the transition 
point -Ffkt as the value of F where the slope of v(F) 
has a maximum. After the transition point the system is 
in a fluid-sliding state [see Fig. |(b)]. All kinks (and an- 
tikinks) have disappeared, and the system is completely 
spatio-temporally chaotic. A short chain with only a few 
kinks still stays in the kink-dominated regime beyond 
Ffkt- Eventually, it jumps to the fluid-sliding state or 
to the solid-sliding state (see Fig. ^). A hysteresis occurs, 
and the chain goes back to the kink-dominated state at 
F ps Ffkt- 




FIG. 10. The velocity distribution in the kink-dominated 
regime near the transition point. The statistics is obtained 
from 5000 samples. Each sample is the sliding velocity v 
averaged over 200 time units. The velocity of the quiescent 
phase is given by vo w 0.1457. The numbered vertical lines 
denote the number N p of kink-antikink pairs. The parameters 
are the same as in Fig. ^. 

At the transition point the sliding velocity strongly 
fluctuates. These fluctuations set in already much below 
the transition point. They lead to a larger value of the 
average sliding velocity compared to the value for small 
chains. A typical example is shown in figure ||. One sees 
bursts of activity above a level given by the value of v for 
small chains (see Fig. |7|). A detailed look onto the dy- 
namics of the chain reveals that an increase of v is caused 
by the production of kink-antikink pairs |lO|Jl4|Jl^| . These 
pairs usually appear behind a 2-7r-kink cluster. This may 
be the reason why for small chains the kink-dominated 
states survive beyond the transition point because the 
probability for a 27r-kink cluster is too small. Each kink 
and antikink contributes to the sliding velocity of the 
chain. That is, v is given by 



„ M + 2N p 
2ir -Cu, 



(10) 



where is the velocity of the kinks and antikinks and 
N p is the number of kink-antikink pairs. In the quies- 
cent state without kink-antikink pairs, the average slid- 
ing velocity is given by vq — 2irMck/N. Thus the ratio 
v/(vq/M) gives directly the total number of kinks and 
antikinks (i.e., M + 2N P ). Figure |i"o| depicts the distribu- 
tion of the sliding velocity shown in Fig. ^. In addition 



to the large peak at the quiescent state one sees clearly 
small peaks for N p = 1 and 2. 

For a approaching zero (or any other integer multiple 
of 2ir) the average sliding velocity of the kink-dominated 
sliding state also approaches zero. It disappears at a — 
0. Thus, the system goes from the fluid-sliding state 
directly to a stationary state when F is decreased below 
the transition point Ffkt- If the value of b is not too 
large [i.e., b — 0(1)], the stationary state will be the 
ground state jl9). The transition is often accompanied 
by transients where a few kink-antikink pairs survive for 
a considerably long time. But in all simulations these 
pairs eventually disappeared. 

We found that the transition point -Ffkt is nearly in- 
dependent of a (as long as a/2ir is near an integer value). 
Table |] shows that -Ffkt increases with b weaker than 
linearly. 

In order to understand the transition between the fluid- 
sliding state and the kink-dominatcd sliding state we 
compare the FK model with a simpler model, namely, 
one particle in a tilted spatially periodic potential plus 
additive white noise. This model was studied in details 
by Risken and VoUmer |^(|. The single particle and 
the noise correspond to the center of mass of the FK 
model and the chaotic motion of the internal degrees of 
freedom, respectively. Of course the noise is neither ad- 
ditive nor white. Its strength depends on the state. It 
is obvious that the solid-sliding state and the stationary 
state of the FK model correspond to the running state 
and the locked state of the simple model in the absence 
of noise. We suggest that the fluid-sliding state and the 
kink-dominated state of the FK model correspond also 
to the running and locked states of the simple model 
but now with noise. Risken and Vollmer showed that the 
bistability between the running state and the locked state 
disappears even for infinitesimally weak noise. There is 
a well-defined transition point F2 which is smeared out 
for increasing noise strength. For 7 <C \fb it is given by 



F 2 = 3.35767 A 



(11) 



Thus we expect Ffkt ~ F2. Table | shows that this is 
indeed the case especially for small values of b. 

TABLE I. The transition point Ffkt between fluid-sliding 
states and kink-dominated sliding states and the transition 
point F2 between the running state and the locked state of sin- 
gle particle in a periodic potential under the influence of weak 
noise. The parameters are N = 500, M = 0, and 7 = 0.05. 
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0.096 
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0.187 
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0.444 
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0.1679 
0.2374 
0.2908 
0.3358 
0.3754 
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III. NONEQUILIBRIUM FREEZING AND 
MELTING 

The aim of this section is to take a closer look at the 
hysteresis loop between stationary states and the fluid- 
sliding state (see Fig. ^). There are two different types 
of transitions. A depinning-pinning transition where the 
(chaotically) sliding chain turns into a stationary one. ft 
is a kind of nonequilibrium freezing and, as in usual first- 
order phase transitions, the chain can be "supercooled" 
below the transition point Fdp ■ That is, for F slightly be- 
low Fdp the chain does not freeze immediately. It takes 
a while until a critical nucleus has appeared. The sec- 
ond transition is the pinning-depinning transition. The 
transition point Fpd depends on the stationary state and 
therefore on the history of the system. In an actual ex- 
periment where one sweeps through the hysteresis loop 
at a finite rate, one will therefore not get well-defined 
transitions points but more or less broad distributions. 
An example of such a (numerical) experiment is shown 
in the inset of Fig. 0(a). It is typical that the distribution 
of Fdp is narrower than the distribution of Fpo- For a 
quasistatic sweep the distribution of Fdp becomes sharp, 
whereas the width of the distribution of Fpd is nearly 
independent on the sweeping rate. 

From equilibrium thermodynamics it is well-known 
that melting and freezing occur at the same tempera- 
ture and the transition is of first-order. In our case far 
from thermal equilibrium the situation is different. In 
the overdamped limit the pinning-depinning transition 
point Fpp) is identical with the depinning-pinning transi- 
tion Fdp but the transition is of second order (|^J2^] . In 
the underdamped case we have bistability between slid- 
ing states and stationary states, i.e., Fdp < Fpd- 
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FIG. 11. Nucleation of a stationary state in the depin- 
ning-pinning transitions. The time step between two snap- 
shots is St — 4n/v. The snapshot just before nucleation is 
highlighted. The parameters are N = 500, M = 160, 6 = 2, 
7 = 0.05, and F = 0.097. 

A further characteristic feature of nonequilibrium 
melting and freezing is the change of the effective tem- 
perature [see e.g. Fig. |l](b)]. This is a general prop- 
erty which is also found in similar models with nonzero 
temperature of the environment [ p!o|l^Jl^ |. That is, the 



fluid-sliding state has an effective temperature which is 
independent on the temperature of the environment as 
long as the latter one is considerably less than the for- 
mer one. The effective temperature is a measure of the 
enhanced energy flow due to the excitation of phonons. 

A. Freezing: The depinning-pinning transition 

We have studied in detail the depinning-pinning tran- 
sition from a fluid-sliding state to a stationary state. We 
have chosen a value of a where the fluid-sliding state is 
completely spatio-temporally chaotic. Like in an ordi- 
nary first-order phase transition at thermal equilibrium 
the transition is caused by a nucleation process. That 
is, a small portion of the chain becomes stationary, and 
the fronts between this stationary nucleus and the slid- 
ing chain propagate into the chain. An example of such 
a nucleation process is shown in figure |ll|. One sees that 
states appear behind the fronts which can be character- 
ized by an average particle distance a (or density 1/a) 
which is roughly constant. Note that the values of a be- 
hind both fronts have to be different. This can be under- 
stood by the following argument. In the previous paper 
we have seen that because of the conservation of the num- 
ber of particles the velocity of a front traveling from one 
particle to the next is given by c = (v\ — i>2)/(a2 — &i), 
where the average particle distance and the average slid- 
ing velocity on both sides of the front are given by ai/% 
and Ui/2, respectively. The chaotic sliding state is char- 
acterized by ai = 2irM/N and v\ > 0. For the stationary 
state V2 = holds. In order to have fronts traveling in 
opposite directions (see Fig. [TT|) the average particle dis- 
tances of the stationary states have to be different (one 
larger than a and one smaller than a). We found that 
they are always the nearest integer multiples of tt, that is, 
for < a < 7r they are zero and n. From this considera- 
tion it is clear that after the depinning-pinning transition 
the system cannot be in the ground state |ljj which is 
characterized by a single domain with a uniform a. In- 
stead the chain will be in a stationary two-domain state 
which is not quite perfect because each domain contains 
defects at low density. 

A similar depinning-pinning transition was found in 
a seemingly completely different model, namely, the 
quenched Kadar-Parisi-Zhang equation with negative 
nonlinearity p3| . It is a nonlinear diffusion equation 
which models the surface growth in disordered media. 
Here the dissipation is strong and the randomness of the 
pinning landscape is important contrary to the weakly 
damped FK model. Nevertheless the pinning-depinning 
transition is of first order and the state which freezes out 
of the sliding state is a two-domain state, as in the FK 
model (compare Fig. 2 of Ref. with Fig. |il| ). 

The duration of freezing t]j is the sum of the nucleation 
time tjv and the growth time to- The nucleation time 
is the time that evolves until a critical nucleus appears. 
A critical nucleus is a nucleus which is large enough to 
grow into the fluid-sliding state. The nucleation time t^ 
will be a Poisson distribution if the probability for the 
appearance of a critical nucleus in a short time-step is 
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small, 



N) 



(12) 



It is well-known that the average and the standard de- 
viation of a Poisson-distributed value are identical, i.e., 
(iiv) = Atjv = 9. For systems much larger than the 
critical nucleus, the probability for nucleation increases 
linearly with the system size, i.e., 9 ~ 1/N. 

The time to for a critical nucleus to grow up to the 
stationary state does not fluctuate as strongly as tjv be- 
cause it is roughly given by the system size divided by 
the sum of front velocities. Thus (to) ~ N. 
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FIG. 12. Cumulative density of the duration times to- 
Each curve is obtained from 500 numerical nucleation experi- 
ments. The parameters are a = 167r/25, b = 2, 7 = 0.05, and 
F = 0.097. 

Figure [l2] shows that the nucleation process is indeed 
characterized by a Poisson distribution with 9 ~ 1/iV. 
From fits of the curves shown in Fig. [l^ we found 9N = 
(2.1±0.1)T0 6 . One can also see a shift of the distribution 
for increasing N to larger values of tjj. This reflects the 
fact that (to) increases with N. 

For values of a between 2.45 and 3.09 (and b = 2, 
7 = 0.05, N = 500, and F « F F kt) we found that the 
fluid-sliding state changes its character. A domain ap- 
pears where the particles are stationary (with two parti- 
cles per potential well). This domain which is surrounded 
by spatio-temporal chaos has a constant size and travels 
through the chain with a constant velocity. The transi- 
tion from this so-called traffic-jam state |Q to a station- 
ary one occurs also via nucleation. The critical nucleus 
seems to appear always at the back of the stationary do- 
main. It reverses the propagation direction of the front. 
That is, the chaotic state stops traveling into the sta- 
tionary one. Instead a front propagates into the chaotic 
state leaving behind a stationary state different from the 
already existing one. Thus the result is again a station- 
ary two-domain state where one domain already exists 
at least partially before the transition. The nucleation 
probability is roughly independent of the system size N 
because the nucleation site is predetermined. 
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FIG. 13. An example of a micro-slip. Snapshots at an 
interval of 10 time units are shown. The inset shows Xj for 
j <E [160, 233]. The parameters are N = 233, M = 89, b = 2, 
7 = 0.05, and F from 0.11554 until 0.115575 at a rate of 10~ 7 . 



B. Melting: Pinning-depinning transition 

The transition from a stationary chain to a sliding one 
is the pinning-depinning transition. The transition point 
Fpd depends on the stationary state. It corresponds to 
a saddle-node bifurcation where the particular station- 
ary state annihilates with an unstable stationary state. 
In the overdamped limit Fpn is uniquely defined by the 
saddle-node bifurcation of the last stable stationary state 
which is usually the state which develops adiabatically 
out of the ground state for F — 0. The disappearances 
of the other stationary states in a saddle-node bifurca- 
tion lead only to more or less local rearrangements of the 
chain. We call such a rearrangements a micro-slip. A 
micro-slip changes the center of mass of the chain but 
does not lead to sliding. An example is shown in fig- 
ure |l^. This behavior and the statistics of micro-slips 
has been studied mainly in models with random external 
potential ||,|||. 

In the underdamped regime the energy gained from a 
local rearrangement is not dissipated immediately. This 
may lead to an avalanche which turns the whole sys- 
tem into a sliding state. Whether a saddle-node bifur- 
cation leads to a micro-slip or to a transition to sliding 
depends on the stationary state and on the damping con- 
stant 7. For each stationary state one can presumably 
find a critical value 7 C which distinguishes both cases: 
For 7 > 7 C we get a micro-slip, otherwise a transition 
to sliding. In our simulations we found only a few (not 
more than three) micro-slips for 7 = 0.05 and b — 2. The 
pinning-depinning transition point depends strongly on 
the history of the system because the actual value of Fpjy 
depends on the stationary state. Sweeping several times 
through the pinning-depinning hysteresis loops, we find 
a broad distribution of Fpd even for very small sweep 
rates. 
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IV. SUMMARY AND CONCLUDING REMARKS 

In this paper we have shown that in the weakly damped 
and strongly driven FK model spatio-temporal chaos ap- 
pears. This chaotic state, called fluid- sliding state, has an 
effective temperature which reflects the fact that phonons 
are excited. This has to be contrasted by the solid-sliding 
state where no phonons are excited pq| . The excitation 
of phonons opens up additional channels of dissipation. 
Thus friction in the fluid-sliding state is larger than in 
the solid-sliding state. 

The velocity-force characteristic shows a pronounced 
transition if the average particle distance a is near but 
not identical to an integer multiple of the period of the 
external potential. It is a transition between the kink- 
dominated sliding state and the fluid-sliding state. Below 
the transition point -Ffkt sliding is caused by the prop- 
agation of 2tt kinks. Approaching the transition point 
from below leads to an increasing production of kink- 
antikink pairs. Above the transition point all kinks and 
antikinks disappear and the chain is in the fluid-sliding 
state with an average sliding velocity which depends only 
weakly on the average particle density 1/a. In the kink- 
dominated regime the average sliding velocity is propor- 
tional to amod27r. The transition point -Ffkt is roughly 
independent of a. For small values of b it is very well 
approximated by the transition point between the run- 
ning and locked solution of an (N — 1)-FK model with 
infinitesimally small additive white noise [i.e. Eq. (Ell)]. 
This raises the question whether the dynamics can be 
reduced to a center-of-mass motion plus some nontrivial 
(e.g. colored, state-dependent) noise term. 

The noncquilibrium freezing (i.e., the transition from 
sliding to stationarity) of the fluid-sliding state is a nu- 
cleation process. The resulting stationary state has two 
domains of different particle densities. 

The nonequilibrium melting (i.e., the transition from 
a stationary state to sliding) of such a stationary state 
corresponds to a saddle-node bifurcation where the sta- 
tionary and, of course, the stable state annihilates with 
its unstable counterpart. In the case of nonzero envi- 
ronmental temperature the transition occurs a bit earlier 
because thermal activation overcomes the barrier which 
decreases to zero at the saddle-node bifurcation. Not 
all saddle-node bifurcations lead to melting. They may 
lead to a local rearrangement of the particle configuration 
called micro-slip. The pinning-depinning transition point 
depends on the history because each stationary state has 
another bifurcation point. 

Our results are qualitatively very similar to results 
of similar models. Persson studied a two-dimensional 
Lennard-Jones liquid on a corrugated potential with 
square symmetry at a nonzero temperature Jl2| . He var- 
ied the temperature between 1/3 and 1/2 of the melting 
temperature (the thermal energy was roughly a tenth 
of activation energy for single particle diffusion). He ob- 
tained velocity-force characteristics and effective temper- 
ature plots similar to Fig |l|. Because of the relatively high 
temperature of the environment he did not find hystere- 
sis between the solid-sliding state and the fluid-sliding 
state. As in the FK model, the transition from the fluid- 



sliding state to a stationary state is also a nucleation 
process. Persson found that it occurs when the effective 
temperature is equal to the melting temperature at ther- 
mal equilibrium. Because the velocity distribution of the 
fluid-sliding state was Gaussian, he argued that the fluid- 
sliding state is in a kind of thermal equilibrium. There- 
fore it has to freeze bel ow th e melting temperature. But 
as we have seen in Sec. II A, Gaussian distributed veloc- 



ities do not imply a quasi-thermal equilibrium. It would 
be interesting to look whether in Persson's model the 
equipartition theorem of the phonon modes is fulfilled or 
not. This is clearly a better test on thermal equilibrium. 

Granato et al. studied a two-dimensional Frcnkcl- 
Kontorova model at nonzero temperatures |ll| . The ex- 
ternal potential is corrugated only in the direction of the 
applied force. Also the particles can move only in this di- 
rection. They did the simulations at a temperature which 
corresponds to 1/4 of the activation energy and which is 
roughly 1/5 of the melting temperature. They also found 
a behavior like in Fig. |l| again without a hysteresis be- 
tween the solid-sliding state and fluid-sliding state. They 
also confirmed the observation of Persson that the tran- 
sition from sliding to stationarity occurs at the melting 
temperature. 

Braun et al. investigated a generalized FK model 
that is quite similar to the model studied by Persson 
0,0Jl§. The main difference is that the substrate po- 
tential is anharmonic and the interaction potential is an 
exponential repulsion. They did the simulations mainly 
for two different temperatures which correspond to 10~ 3 
and 1/20 of the activation energy for single particle dif- 
fusion. They report results for one-dimensional chains of 
atoms as well as for two-dimensional layers. There seems 
to be no qualitative differences between one and two di- 
mension. Again velocity-force characteristics and effec- 
tive temperature plots are similar to our results. They 
found a hysteresis between solid sliding and fluid slid- 
ing. Its width decreases with the environmental temper- 
ature and disappears eventually (at a temperature which 
is roughly 1/5 of the activation barrier for hopping of 
uncoupled single particles). Near fully commensurate 
particle densities they found a transition (denoted by 
-Fpair) which is similar to the transition from the kink- 
dominated sliding regime to the fluid-sliding regime of 
the FK model. Because their chain was short (i.e., N = 
105) this transition shows hysteresis for very low tem- 
peratures. Instead of a fully chaotic fluid-sliding regime 
they always found two-domain states with alternatively 
running and locked particles (traffic-jam regime). They 
also measured the transition point -F pa ir as function of 
the damping constant p8[ . It increases roughly linear 
with the damping constant but the authors seemed not 
to be aware that ([H]) is again a remarkable good approx- 
imation (i.e. errors are less than 10%) for F pa i r . 
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